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Abstract. We present a new proof of a theorem of Mallat which describes a con- 
struction of wavelets starting from a quadrature mirror filter. Our main innovation 
is to show how the scaling function associated to the filter can be used to identify 
a certain direct limit of Hilbert spaces with (R) in such a way that one can im- 
mediately identify the wavelet basis. Our arguments also use a pair of isometrics 
introduced by Bratteli and Jorgensen, and exploit the geometry inherent in the 
Cuntz relations satisfied by these isometrics. 

A wavelet is a function i) e such that 

{V^j- fc : X ^ 2-^l^^{2-^x -k): i,ke Z} 

is an orthonormal basis for There are, remarkably, many different wavelets, 

and they have proved to be enormously useful in both theory and applications. So 
there has been a great deal of interest in methods of constructing wavelets. One 
famous construction of Mallat [Zj starts from a quadrature mirror filter: a function 
mo : T — >• C such that 

|mo(-2)P + |mo(— = 1 for every z eT. 

Our goal here is to present a new proof of Mallat's theorem based on the concept of 
a direct limit. 

Mallat proved his theorem in two stages. From the filter mo he built a multires- 
olution analysis, in which a central role is played by a scaling function (p G L^(M) 
satisfying (j){2x) = mo(e^'^"')(/)(a;) [71 Theorem 2]. He then used what he described as a 
"by now classic" algorithm to generate the wavelet [3 §4]. Mallat's construction has 
since been refined and discussed in several books. For example, ^ §5.3-4] contains 
a relatively elementary proof of his theorem, in which some of the analysis has been 
simplified but the overall strategy is that of Mallat. In our proof, the scaling function 
still plays a central role: we use it to identify a certain direct limit with L^(M), and 
the existence of the wavelet then follows almost immediately from the geometry im- 
plicit in some operator-theoretic equations called the Cuntz relations. The analytic 
content of our proof is much the same as that in the standard sources, and we refer 
to them for details, but our organisation seems to be quite different. 
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Our arguments seem to be more natural in the Fourier or frequency domain, so 
we work there throughout, and our construction yields the Fourier transform of the 
wavelet. One effect of working in the Fourier domain is that the scaling equation 
(Equation below) involves multiplication rather than convolution. In the interests 
of clarity, we shall consider only the classical (dyadic) wavelets, but we are optimistic 
that our approach will also shed light in other situations where wavelet bases are 
used. 

Throughout, mo will be a quadrature mirror filter such that mo is smooth at 1, 
mo(l) = 1 and mo{z) ^ for z in the right half-circle. We define mi : T ^ C by 



rax{z) := zttiq^—z). 
We now define two operators Sq and on L^(T) by 

(1) {SJ)iz):=2'/'m,iz)fiz'). 

Our starting point is the following observation of Bratteli and Jorgensen P] : 
Lemma 1. The operators Si satisfy SqSq = 1 = S^Si and SqSq + SiSl = 1. 
To prove this, first verify that the adjoints S* are given by 

(5;/)(e2"") = 2-1/2 ^ ^^^^(^/(e""0 +^^i^Fe^/(-e""')), 

and then compute S*Si and SqSq + SiSl- 

The formal computation in Lemma ^ has some very interesting geometric conse- 
quences. The relations S*Si = 1 say that the operators Si are isometrics of L'^{T) 
into itself, and imply that the operators SiS* are the orthogonal projections onto the 
ranges of the Si (which are automatically closed because the Si are isometrics) . Since 
the sum of two projections is a projection only when their ranges are orthogonal, 
the Cuntz relation SqSq + SiSl = 1 implies that the ranges SiL^(T) of the Si are 
orthogonal complements of each other, so that 

(2) L\T) = SoL'^iT) © SiL\T). 

We now recall a construction from algebra. Suppose that we have Hilbert spaces 
Hn and isometrics T„ : Hn Hn+i for all n G N. The direct limit {Hoc, Un) consists 
of a Hilbert space Hoc and isometrics f/„ : Hn — ^ Hoo which satisfy Un+i o Tn = Un 
and which have the following universal property: for every family of isometrics {Rn} 
of Hn into a Hilbert space K such that Rn+i oTn = Rn, there is a unique isometry 
Roo '■ Hao — > K such that R^o o Tn = Rn for every n. We illustrate this universal 
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property with the diagram: 



Uo 




The uniqueness imphes that H^o = U^i UnHn, since otherwise we could define i?oo 
arbitrarily on ((J UnHn)'^. 

It is not hard to see that the direct limit exists. Indeed, identifying each h G 
with all its images Tm-iTm-2 ■ ■ ■ Tnh G Hm defines an equivalence relation ~ on 
the disjoint union |J Hn, and because the T„ are isometrics and hence inner-product 
preserving, the set H' :— (|J Hn) / ~ of equivalence classes is naturally an inner- 
product space; completing H' gives a Hilbert space H^o, and the maps f/„ which send 
elements of Hn to their class in H' C H^o have the required properties. However, 
the important point is that the construction does not matter, since the universal 
property identifies the direct limit up to isomorphism: to identify with a Hilbert 
space K, for example, we just need to find isometrics Rn : — > X as above such 
that K = RnHn, and then R^o is an isomorphism of ifoo onto K. 

If we start with a single isometry S* on a Hilbert space H, we can take the direct 
limit (i^oo, Un) of the system in which every Hn is H and every T„ is S. Now consider 
the diagram 




where the horizontal arrows going into Hf^, are there to remind us that we have 
isometrics Un of every copy Hn of H into i?oo- Applying the universal property of 
the top row to the downward arrows (or more properly, to i?„ := ?7„ o S) gives an 
isometry 6*00 on H^o which is characterised by SooiUnh) = Un{Sh). Similarly, applying 
the universal property of the bottom row to the NE arrows gives an isometry Iqo which 
is characterised by loo{Unh) = Un+ih. Then for every n and every h & H we have 

{looSoc){Unh) = loo (5*00 t/„/i) = ioo{Un{Sh)) = Un+i{Sh) = C/^/i and 
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which imphes that loo is an inverse for Soo- This process of passing to the direct hmit, 
therefore, turns the isometry S into a unitary S^o- Since 

(3) Soo{Un+lh) = Un+l{Sh) = UnK 

this unitary is an isomorphism of the copy Un+iH of Hn+i = H onto the copy UnH 
of Hn = H. 

Applying the process described in the previous paragraph to the isometry Sq on 
L^(T) defined in (d)) gives a Hilbert space H^o and a unitary operator 6*00 on H^o- Our 
next task is to identify the direct hmit with L^(M). MaUat proved that, under 
our hypotheses on mo, there is a scaling function} G L^(R) of norm 1 such that 

(4) 0(2x) = mo(e2"")0(x) and 

(5) 5^|0(x + A;)p = l 

for every x G M. Indeed, he proved that the infinite product 



[X) = 

n=l 



JJmo(exp(27ri2-"x)) 



has the required properties (see pages 76-77] or |Hl §5.4]). For example, if mo 
is the characteristic function of the right half-circle, then = X[-i/2,i/2] is a scaling 
function. (In '6, page 225] this is proved under the milder hypothesis that mo satisfies 
Cohen's condition^ which is a necessary and sufficient condition for the existence of 
such a function 0.) 

We now fix n G Z, and define i?„ : = L^(T) L'^iR) by 

iRnf)ix) = 2-"/V(exp(27rz2-"x))0(2-'^x). 

Each Rn is an isometry: indeed, a change of variables and an application of the 
monotone convergence theorem shows that 



\Rnf\\lm) = I 2-"|/(exp(27rz2-"x))0(2-"x)|^rfx 
|/(e2™)0(a:)|2rfx 

= f\f{e''''')\'{j2\^(t + n)\')dt, 

which by (0) is precisely the norm of / in L^(T). The identity (j3)) implies that the 
isometrics Rn are compatible with the isometrics Sq : Hn = -L^(T) Hn+i = L'^{T) 



^In the literature, it is the inverse Fourier transform cj) of this function (/) which is usually called 
a scaling function. 
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in the direct system: 

^n+i(5o/)(a;) = 2-("+')/'2i/2rno(exp(27ri2-('^+^)a;))/(exp(27ri2-("+i)a;)2)</.(2-("+i)x) 
= 2-"/V(exp(47rz2-("+i)x))0(2(2^("+^)a;)) 

= iRnf)ix). 

Thus the universal property of the direct hmit gives an isometry : L^(M). 
Notice that if we set Vn ■= range i?„ = Roo{UnL'^{T)), then the last calculation shows 

that Vn C Vn+l- 

Lemma 2. The isometry R^o intertwines the unitary Soo on and the dilation 
operator D on ^^(R) defined by {D^)ix) = 2^/^^{2x). 

Proof. We let / G = L'^iT), and compute: 

D{R^{Un+imx) = D{Rr,+Mx) 

= 2'/'{Rn+J){2x) 

_2i/22-("+i)/2j(exp(27r22-("+^)2x))0(2-("+^)2x) 

= {Rnf){x) 

= Roo{Unf){x). 

Since we know from Q that o Un+i = Un, this implies that D o R^ and i?oo ° •S'oo 
agree on the range of every f/„+i, and hence on Hoc- □ 

Since 5*00 is an isomorphism of Un+iH onto UnH, it follows from Lemma El that^ 
DVn+i = Vn for every n G Z. The following lemma is proved in [3 pages 78-79] 
and Lemmas 5.47 and 5.48]. (A different argument which proves the analogous 
property of the inverse Fourier transforms is given in Propositions 5.3.1 and 5.3.2 
of 0.) 

Lemma 3. We have Hngz = W ^'^'^ Unez^« = L^(M). 

We write Wn for the complement Vn+i Vn of Vn in Vn+i- Then the subspaces 
Wn are mutually orthogonal, and it follows from Lemma El that L^(]R) decomposes as 
the direct sum ©^ggVTn- Since D^^Vn^i = Vn for every n, we have DWn+i = Wn 
for every n. Thus to find an orthonormal basis for L^(]R), it suffices to find an 
orthonormal basis for one Wn, and then this together with all its dilates will be an 
orthonormal basis for L^(R). 

So we seek an orthonormal basis for Wq. With H = L^(T), we have 

Wo = VieVo = RiH e RqH = r^h e Ri{SoH) = Ri{He SqH). 

At this point we recall from Q that the complement of SqH is the range SiH of the 
other isometry Si. Since RiSi is an isometry, it maps the usual orthonormal basis 

^This formula looks slightly different from the usual dilation property of a multiresolution analysis 
because we are working in the frequency domain; in and 0, for example, the space denoted by 
Vn is the inverse Fourier transform of our Vn ■ 
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{ck : z ^ : k & Ij] for L^(T) into an orthonormal basis for Wq. We deduce that 
the functions 

form an orthonormal basis for Wq. We set 

^{x) := mi(e"")0(2"ix), 
so that the basis elements take the form 

If we now define 

^^.^(x) = {D^iljk){x) = 2^/'^exp{-2TTik2^x)tlj{2^x), 

then {V'i.fc '■ j,k & Z} is an orthonormal basis for L^(]R). 

Since the inverse Fourier transform intertwines D and D^^, and intertwines multi- 
plication by e^'^*'^^ and the translation operator which takes ^ to ^(■ + A;), the functions 

iPj^kix) = 2-^/^{2-^x - k) 

also form an orthonormal basis for L^(M). In other words, is a wavelet, and we 
have proved Mallat's theorem: 

Theorem 4. Suppose that mo is a quadrature mirror filter such that mo is smooth 
at 1, mo(l) = 1 and mo{z) ^ for z in the right half-circle, and let (f) he a function 
satisfying the scaling conditions (jH) and (0) . Define : M ^ C by 

^(x) = e'^*^ mo(-e^^^)0(x/2). 

Then the inverse Fourier transform tp is a wavelet. 
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